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Abstract. A method of radio frequency interference (RFI) suppression in radio astronomy spectral observations 
is described based on the analysis of the probability distribution of an instantaneous spectrum. This method 
allows the separation of the gaussian component due to the natural radio source and the non-gaussian RFI signal. 
Examples are presented in the form of 
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1. Introduction 

The increasing radio frequency interference (RFI) from 
transmissions and spurious signals of commercial telecom- 
munication applications limit the real sensitivity of radio 
telescopes at certain frequencies. This radio-ecological en- 
vironment at radio observatories appears to be worsening 
each year, while the observing systems get more sensi- 
tive. Further improvements in antenna quality and in re- 
ceiver parameters do not give the expected results because 
of limitations due to RFI. Therefore, it is very difficult 
to reach the design sensitivity, especially at low frequen- 
cies, even for radio telescopes situated in secluded areas 
(SpoelstrCEMI KahlmannQIMl Cohen [TH|) . 

Several methods of RFI excision have been proposed 
during recent years. The successful application of these 
methods depends on the type of radio telescope, type of 
observations, type of RFI, and there is no universal RFI 
mitigation algorithm. 

1. Type of radio telescope - single dish, connected 
radio intererometer, VLBI. Single dish radio telescopes are 
especially vulnerable to RFI. Radio interferometric sys- 
tems suffer from RFI in a less degree, (Thompson 1982), 
but auxiliary noise power measurements made at each an- 
tenna for calibration purposes are distorted by RFI, as for 
the single dish case. 

2. Mode of observations - continuum or spectral. 
In continuum observations it is possible to sacrifice some 
part (5-10%) of the data stream polluted by a strong RFI 
and save the remained part without significant loss of sen- 
sitivity. In spectral observations RFI and signal-of-interest 
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may occupy the same spectral region, so it is impossible 
to delete this particular part of the spectrum. 

3. Type of RFI - impulse-like bursts, narrow-band 
or wide-band. Theoretical consideration and experimen- 
tal data show that RFI may be basically represented as 
the superposition of two types of waveforms: impulse-like 
bursts and long narrowband and strongly correlated ran- 
dom oscillations (Middleton 119721 113771 Lemmon [TW|) . 
RFI of the first type can be effectively suppressed using 
fast real time thresholding before averaging at the total 
power detector or correlator (Fridman 119961 Weber et al. 
1997). This is the processing in temporal domain. 

RFI of the second type should be processed in the fre- 
quency domain, where there is a good contrast between 
a quasi- continuum spectrum of the system noise plus ra- 
dio source signal and the narrow-band RFI signal. RFI in 
continuum observations may be excised using adaptive fil- 
tering methods (Fridman 1998), such that the parts of the 
spectrum affected by RFI are used with a smaller weight- 
ing than the other "unpolluted" spectral regions during 
the averaging in the time-frequency plane. This sort of 
processing cannot be applied in spectral line observations, 
where RFI may contaminate the same spectral region that 
is of interest for radio astronomer. In (Ellington, Bunton 
and Bell EOOO]) RFI from navigational satellites GLONASS 
was parametrically modelled and the model of RFI signal 
was subtracted from the radio telescope output so that 
the 1612MHz OH spectral line was not distorted. 

Adaptive interference cancelation (similar to adaptive 
noise cancellation, Widrow &Stearns 1985) can be use- 
ful for single dish observations with an auxuliary refer- 
ence antenna, pointed "off-source". The use of an auxil- 
iary reference channel provides some benefits (Barnbaum 
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& Bradle\ fl998p . but not all radio telescopes are or can 
be equipped with such a reference channel. Furthermore, 
the effectiveness of this adaptive noise canceling strongly 
depends on the sensitivity of the reference channel. The 
reference channel should at least be equally sensitive than 
the main channel for successful RFI suppression (Fridman 
11557112000) . 

A new method of post-correlation RFI removing was 
proposed in (Briggs, Bell and Kesteven [2000). Cross- 
correlation products of four signals (two main channels 
and two reference channels) are processed, using closure 
complex amplitude relations. RFI signal is estimated and 
subtracted from the outputs of main channels. So RFI sup- 
pression can be made off-line with the averaged data, see 
also (The Elizabeth and Frederick White Conference on 
Radio Frequency Interference Mitigation Strategies [1999]). 

Methods of spatial filtering in the temporal and fre- 
quency domains applicable for radio interferometers are 
discussed elsewhere (Leshem & van der Veen [T999). 
Spatial filtering methods use difference in the direction-of- 
arrival of radio source and RFI. The RFI emission of well 
spatially localized sources could be suppressed in multiele- 
ment radio interferometers using adaptive array philoso- 
phy, when zeros of a synthesized antenna pattern coincide 
with directions-of-arrival of undesirable signals (adaptive 
nulling techniques). This approach is popular while con- 
sidering new projects of superlarge phased array radiote- 
lescopes, (van Ardenne 119961 Bregman I2000|) . However, 
there are some limitations in using this techniques for 
the existing large radio interferometers like WSRT, VLA, 
GMRT: very sparse antenna arrays, only-phase computer 
control, complex impact on image postprocessing. Spatial 
filtering is effective when RFIs are strongly correlated at a 
radio telescope antennas (like in a half-wavelength phased 
array). RFI at the adjacent antennas (hundreds meters 
separation) are highly correlated, but when the distance 
between antennas is more than 1km (and especially for 
bases 10-100km) RFI at these points are less correlated, 
and effectiveness of the adaptive nulling is not so high, 
as for the continuous-like half-wavelength phased arrays. 
RFI at the sites of VLBI system, separated by hundreds 
and thousends of kilometers, are practically uncorrelated. 
So the impact of an RFI is an increased variance at corre- 
lator output, which is noticeable in the case of strong RFI 
(comparable with the system noise power). 

In this paper, I consider a method of RFI excision 
based on the analysis of the probability distribution of the 
power spectrum at the radio telescope output. System and 
radio source noise generally has a gaussian distribution 
with a zero mean. If a Fourier transformation is applied 
to such an "ideal" signal, the real and imaginary compo- 
nents in every spectral bin are also gaussian random val- 
ues with zero mean. The instantaneous power spectrum, 
which is a square of magnitude of the complex spectrum, 
has an exponential distribution (or, in other terms, a chi- 
squared distribution with two degrees of freedom). Any 
changes of the input signal of this model, due to the pres- 
ence of an RFI, yield a change of its statistics. The power 



spectrum in this case has a non-central chi-squared dis- 
tribution with two degrees of freedom. Radio astronomy 
observations generally employ the averaging of the power 
spectrum over some integration period, which results in a 
conversion of the probability distribution of the averaged 
power spectrum into a gaussian distribution. 

Analysis of sample probability distribution before av- 
eraging allows an effective separation of these signals by 
means of real-time processing facilities employing digital 
signal processing. The method may be effective for single 
dish observations without reference channel. The descrip- 
tion of this separation by using the analysis of the power 
spectrum distribution will be presented below. I will also 
present some computer simulation examples, and some re- 
sults from test observations at the Westerbork Synthesis 
Radio Telescope (WSRT). 

2. Power spectrum probability distributions and 
moments 

A radio telescope receiver output signal is described as 
follows: 

r(t) = x sig {t) + x sys (t) + x RFI (t), (1) 

where x sys (t) is a system noise signal (receiver, feed, an- 
tenna), x s ig{f) is the radio source noise signal with some 
spectral features to be estimated, and XRFi{t) is the 
RFI waveform. x sys (t) and x S i g {t) are each supposed to 
be random processes with gaussian probability distribu- 
tions, zero mean and variances <j 2 sys and cr 2 ig , respec- 
tively. However, XRFi(t) is a narrow-band interference sig- 
nal with variance o RFI . After Fourier transformation on 
a time interval T, we get the instantaneous complex spec- 
trum from eq.(l): 

S(f)= [ r(t) cos 2n f tdt + j [ r(t) sin 2nftdt = 
Jo Jo 

= S.re(f)+jSJm{f) (2) 

Probability distributions of the real part of the spec- 
trum Sjre(f) and the imaginary part SJ,m(f) are also 
gaussian with a mean A(f) and a variance <r^(/) + A(f) 2 l , 
where component A(f) is determined by the RFI and the 
cr 2 (f) is determined by system + signal spectral densities. 
A power spectrum (periodogram) is written as: 

P(/) T = S(f)S(f), (3) 

which corresponds to an envelope quadratic detector op- 
eration, the overbar denotes complex conjugation. Here 
P(f)r at each frequency bin / has the probability distri- 
bution of the square of the envelope of a randomly phased 
sine wave and a narrow-band gaussian process (noncentral 
X 2 distribution with two degrees of freedom): 

1 P + A RFI A„^t^/~P 

*(P) = i-^I„«,P> 0, (4) 
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Fig. 1. The non-central % 2 probability distribution as a 
function of a)(A RF j /a n ) and b) the order M. 

which in the absence of RFI is an exponential distribution: 



(5) 



After averaging M periodograms (3), P a ver(f) 



I2 m =i p ™(f), we have for P a ver the non- central \ dis- 
tribution with 2M degrees of freedom, fWhalen ll971[) : 



W(P )- 1 f P °n 



2a 2 ■ 



IM 



.A RFI y/PM 
-H ~2 ) 



and in the absence of RFI: 



Wo(Paver) 



(2^) M r(M) 



pM-l e -P 2 /2al 



(6) 



(7) 



When the noncentral parameter A RF j/a n > 1, both 
distributions (4) and (6) tend to acquire a gaussian form. 
For big M (central limit theorem) the distribution (6) 
also is transformed into a gaussian. This near-gaussian 
distribution characterizes the power spectra at the out- 
puts of common spectrum analyzers. Figure la shows how 
dramatically the distribution (4) is transformed, while 
Ab.fi/ o~ n grows from zero to 5. Figure lb represents 
this kind of transformation with periodogram averaging 
applied. Figure 2a shows several sections of the three- 
dimensional presentation of Fig. la. 

One observes that the presence of RFI changes the 
form of the distribution w(P) of the unaveraged peri- 
odogram. This property may be used for separating the 
(system + signal) noise component <r^(/) from the RFI 
component A RF i{f). 

The most promising method is to obtain an estimate of 
the sample probability distribution w[P(f)] (histogram), 
and subsequently estimate the S^(/) and A RF j(f), which 
best fit this sample distribution (symbol^denotes sample 
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Fig. 2. a) The non-central \ 2 probability distribution as 
a function of a) the value of the noncentral parameter 
(Arfi/o-) and b)the skewness and excess as a functions 
of A RFI /a„ . 

average). However, performing this processing in real-time 
with a sufficiently wide band (at the MegaHertz level) is 
a difficult task even for the most powerful modern Digital 
Signal Processors (DSPs). 

A more attractive practical approach is to use the sam- 
ple moments of the power spectrum. The first four statis- 
tical parameters, mean, variance, skewness and excess, for 
distribution (4) are: 



skew = 



mean = 2er* 
var — 4cr^ 
2al 



A 2 



A 2 rfA 3/2 



6(< 



^^RF^TJT 



AIfi) 2 ' 



(8) 
(9) 

(10) 
(11) 



In the absence of RFI, the mean=2a 2 , the variance 
= 4<x^, the skewness = 2, and the excess = 6. For strong 
RFI, the skewness — > and the excess — ¥ 0, which is 
a good indicator of the presence of RFI. Fig. 2b shows 
this dependence of skewness and excess on the parameter 
Arfi /On- 
Having the sample values of mean, var, skew, exc, this 
redundant system of nonlinear equations can be solved 
with respect to o~\ and A 2 RFI , with the obvios constraints 
a 2 > 0,A RFI > 0. Analysis of the "solution surface" 
shows that the solution is unique. Even on the basis of 
only the first two equations, it is easy to derive that at 
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Fig. 3. A block diagram of the computer simulations. 
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The value 2a 2 i (f) as a function of frequency provides 
"clean" power spectrum. 

Sometimes radio astronomers use a cross-correlation 
spectrum instead of the auto-correlation spectrum, with 
the aim to suppress large system noise spectral compo- 
nent. The probability distribution and the statistical pa- 
rameters for this case are given in the Appendix. 

3. Computer simulations 

Some computer simulations have been done with the aim 
to investigate how effective the separation of the gaussian 
noise component and the non-gaussian RFI can be. Figure 
3 represents a block diagram of the simulation program. 
The three signal components are as follows: (1) the sys- 
tem noise is a sequence of gaussian random values passed 
through a passband Butterworth filter with cutoff frequen- 
cies [0.1, 0.9] and with order 7, (2) the signal noise imi- 
tating a spectral line is also represented as a sequence of 
gaussian random values passed through a Butterworth fil- 
ter with cutoff frequencies [0.3 0.35] and with order 5, and 
(3) the RFI waveforms consists of a superposition of two 
strong chirp (frequency modulated sinusoidal wave) sig- 
nals. One of the RFI signals coincides in frequency with 
the position of the " spectral line" . These three signals are 
added and applied to a spectrum analyzer (FFT block), 
where an instantaneous complex spectrum is calculated. 
This operation imitates a discrete digital spectral analysis 
with L frequency channels. The generation of the mixture 
noise + RFI and of the Fourier transform was repeated M 
times. The upper part of the block diagram after the FFT 
block corresponds to the proposed RFI excision method, 
where statistical parameters (mean, variance, skewness, 
and excess) are calculated at each frequency bin for the 



M samples. Subsequently the estimates a^(fi), for each 
frequency channel I — 1...L were calculated. 

The lower part of the block diagram in Fig. 3 shows 
the common way of doing spectral analysis: calculating of 
the power spectrum P(f)r = S(f)S(f) and subsequent 
averaging. At this point the two versions of the power 
spectrum <?„(/;) are available: one without RFI excision, 
and another with RFI excision. Figures 4-7 show the re- 
sults of the some computer simulations. Fig. 4a is an ex- 
ample of an input power spectrum without RFI with its 
input signal waveform in Fig. 4b. Fig. 4c shows the "dirty" 
power spectrum without any RFI excision on a logarith- 
mic scale with the level of RFI about 30 dB above the 
system noise. Fig. 4d is the waveform of the input signal 
with a strong RFI source on the same scale as the system 
noise. Finally, Fig. 4e shows the "clean" power spectrum, 
where the "spectral line" is clearly visible and the RFI is 
significantly suppressed, and Fig. 4g shows the "residu- 
als" or signal-of-interest distortions, that is the difference 
between "noRFI" power spectrum a) and "clean" power 
spectrum e). The rms level of these variations is approxi- 
mately at the level of system noise rms ~ a1 ys /^fM in the 
absence of RFI. Fig. 5 gives the spectra of four statistical 
parameters: mean and variance in logarithmic scale, and 
skew ness and excess in linear scale. Fig. 6 and 7 repre- 
sent similar computer simulation, but with the RFI not 
coinciding with the "spectral line". Fig. 7c and d show 
the remarkable characteristics of the skewness and the 
excess: they are "transparent" for the gaussian "spectral 
line" component and react only to RFI presence, that is 
in the absence of RFI, the skewness = 2, and the excess 
= 6, and for strong RFI both skewness and excess — > 0. 
This is in harmony with (10), (11) and Fig. la, Fig. 2b, 
due to "gaussinization" of probability distribution func- 
tion: skewness and excess are equal to zero for gaussian 
distribution. 



4. Test observations 

Test observations were made with the WSRT to demon- 
strate the effectiveness of the proposed method with real 
radio astronomy signals. Digital spectral analysis was per- 
formed with the help of a SIGNATEC PMP-8 board 
bearing eight high performance Texas Instruments DSP 
TMS320C6201 processors (Fridman [50001) • Th e baseband 
(0-1.25 MHz) signal at the analogue output of WSRT fre- 
quency conversion subsystem corresponding to one of the 
radio telescopes (RT6) is digitized in a 12-bit analogue- 
digital converter and then applied to the PMP-8, where 
the FFT and the averaged statistical parameters were cal- 
culated. The final presentation of the "dirty" and the 
"clean" spectra, and of their difference were done on a 
workstation. 

Figures 8, 9, and 10 show the results of three test ob- 
servations of Galactic spectral lines: 21cm HI lines in the 
direction of Crab and Orion, and an 18cm OH line in the 
W3(OH) radio source. An artificial continuous wave signal 
was generated with a Rohdc&Schwarz SMR20 to serve as 
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Fig. 4. Computer simulations: a) the power spectrum with 
spectral line with no RFI; b)the input signal with no RFI; 

c) a two feature RFI power spectrum (logarithmic scale); 

d) thc resulting input signal with RFI; e)the "cleaned" 
power spectrum; g)the difference between a) and e). 




Fig. 5. 

in Fig. 

excess. 



Spectra of the statistical parameters for the model 
4c (with RFI): a)mean; b)variance; c)skewness; d) 



Fig. 6. Computer simulations: a) the power spectrum 
with spectral line with no RFI; b)the input signal with 
no RFI; c) a two feature RFI power spectrum (logarith- 
mic scale); d)the resulting input signal with RFI; e)the 
"cleaned" power spectrum; g)the difference between a) 
and e). 



the RFI signal and was emitted in the direction of RT6 
from the WSRT control building. The figures display the 
data with and without RFI excision applied. These three 
examples represent relatively strong interfering signals in 
astronomical data. The RFI suppression obtained in these 
observational spectra are respectively I7dB in Fig. 8, I3dB 
in Fig.9, and 20dB in Fig.IO. 



5. Conclusions 

1. Single dish spectral observations are most susceptible 
to RFI. Digital spectral analysis with real-time analy- 
sis of the power spectrum higher order statistics pro- 
vides a promising way for narrow-band RFI excision 
taking into account the performance of modern DSP 
systems. The method can be also applied in continuum 
observations. There is no principal bandwidth limi- 
tations, algorithm works in frequency domain, after 
Fourier transformation. Several RFIs at different fre- 
quencies are processed in a parallel way, independently 
of each other, Fig. 4. 
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Fig. 7. Spectra of the statistical parameters for the model 
in Fig. 6c (with RFI): a)mean; b)variance; c)skewness; d) 
excess. 



2. Computer simulations and test observations demon- 
strate the practical effectiveness of this method with 
respect to the narrow-band non-gaussian RFI signals. 

3. The proposed method is a nonlinear procedure and the 
post-processing suppression gain (the RFI suppression 
1 N Rinput 1 1 N Routput ratio) depends on the input in- 
terference/noise ratio IN Ri nput = A 2 RFI /a n in a com- 
plex way. For large INRi nput this gain can be esti- 
mated approximately as 10log(I N Ri nput ) + 5log(M). 
The stronger RFI, the bigger the gain. The results from 
our tests with astronomical data show that the gain 
may range from 13 - 20 dB. 

4. The effectiveness of the method depends also on the 
temporal stability of INR: equations (8)- (11) must be 
solved under the assumption of constant A during the 
averaging. This limits the value for M and hence the 
averaging procedure should be divided on three stages: 
preliminary averaging - RFI subtraction - final averag- 
ing. 

5. The proposed method was tested for the autospectral, 
periodogram analysis, that is the power spectrum and 
its high-order statistics were calculated in real time 
after Fourier transform of the baseband signals. Most 
of the modern radio astronomy spectrum analyzers use 
auto- or cross-correlation techniques with 1-2 bit quan- 
tization and calculate only the first moment of power 
spectrum. I think that the direct periodogram spec- 
tral analysis with 12-16 bit sampling is more robust to 
strong RFI - less harmonic distortion products. This 
is also more suitable for implementation of other RFI 
mitigation methods using the same DSP hardware, as 
it is practiced in the RFI mitigation test system at 
WSRT. 



6. Appendix 

For a product of two gaussian values z = xy with vari- 
ances C7i,cr2, zero means and correlation coefficient p the 
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Fig. 8. RFI excision using the WSRT with bandwidth = 
1.25MHz, sampling frequency = 1.563MHz, and frequency 
resolution = 3.052kHz: a) a neutral hydrogen spectral line 
(21cm) in the Galaxy towards the Crab nebula with RFI 
superposed at 1420.397MHz; b) the same spectrum after 
RFI excision simultaneously made with frame a; c) a mag- 
nification of the spectral line region (channels 30-100) with 
spectra of a) and b) superposed on a record without any 
RFI. 



probability distribution function (PDF) is 



p(z 



--Ko[ 



7TCTlCT2V / ( 1 - P 2 ) ^^(l-p 2 ) 

pz 



exp[ 



C7lCT 2 (l - p 2 ) 

and the corresponding characteristic function is 
f(0 = [1 - 2i£y/wi + e 2 a 1( r 2 (l - p 2 )]- 1 ' 2 . 



(14) 



(15) 



Fig. 11a shows a 3D presentation of how this PDF 
changes with rising p. Fig. lib represents several sections 
of this dependence. It is clear that for p = 1 this distri- 
bution is transformed into the distribution of a squared 
gaussian value. 

For a product of two gaussian values z — xy with non- 
zero means the formula for the PDF is too complex, but 
the characteristic function is simple enough to calculate 
the necessary central moments and hence the skewness 
and the excess. The central moments are more convenient 
to calculate using cumulants: n r = i~ r ^F^/(£)|£=o- For 
variances equal to 1, the means ai, 02, and the correlation 
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Fig. 9. RFI excision using the WSRT with bandwidth = 
1.25MHz, sampling frequency = 1.563MHz, and frequency 
resolution = 3.052kHz: a) a neutral hydrogen spectral line 
(21cm) in the Galaxy towards the Orion nebula with RFI 
superposed at 1420.350MHz; b) the same spectrum after 
RFI excision simultaneously made with frame a; c) a mag- 
nification of the spectral line region (channels 60-120) with 
spectra of a) and b) superposed on a record without any 
RFI. 
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Fig. 10. RFI excision using the WSRT with bandwidth = 
2.5MHz, sampling frequency = 3.125MHz, and frequency 
resolution = 6.104kHz: a) a hydroxyl (OH) emission line 
(18cm) in the Galaxy towards W3(OH) at 1665.8MHz 
with RFI superposed; b) the same spectrum after RFI 
excision simultaneously made with frame a; c) a magnifi- 
cation of the spectral line region (channels 250-300) with 
spectra of a) and b) superposed on a record without any 
RFI. 



coefficient p, the characteristic function is (Deustcbl962]): For p — > 1, we find: 

2iia x a 2 - £ 2 [a 2 - 2pa x a 2 + a 2 ] 
ex Pt ^ r,-, I -PTR oTi • ( 16 ) 



ml cor = A 2 + 2a 2 , 
var cor = 4A 2 <T 2 + 4er 4 , 



(21) 
(22) 



For the cross-spectral density at a particular frequency, we 
have to use /x;2(£) = /(C) 2 due to two degrees of freedom, 
as for the auto-correlation spectrum. Thus the first two 
central moments corresponding to /s 2 (£) and to variances 
a 2 , and a 2 (so that x = Oici, and y = a 2 a 2 ) are: 



mlr 



2pa 1 a 2 + 2xy, 



var cor — 2x 2 o\o 2 + 2y 2 a\u 2 + Apxy<7\<j 2 + 

2a\a\ + 2p 2 a\a\. 



(17) 



(18) 



For <7i — a 2 — cr, and x = y = A/^/2, we have the 
system of equations with respect to a and A: 

ml cor = 2p<7 2 + A 2 , (19) 
Picor = var cor = 2A 2 a 2 + 2pA 2 a 2 + 2a 4 + 2p 2 a 4 (20) 



which is the same as for the auto-correlation spectrum. 
The following formulae are for the third and the fourth 
central moments of the sum of two products: 

Hscor = 4p 3 cr 6 + 12pcr 6 + 6A 2 p 2 a 4 + \2pA 2 (j 4 + 6A 2 a 4 , (23) 

p icor = 24p 4 a 8 + 96pV + 24 2 cr 8 + 48^V^ 6 + 96A 2 p 2 a 6 + 

+ 96 A 2 pa e + 48A 2 cr 6 + 12A 4 p 2 a 4 + 2AA 4 pa 4 + \2A 4 a 4 .{2A) 

Using these formulae for p j , it is possible to calculate skew- 
ness ske 



J cor — 3/2 — 



a 4 V2(2p 3 (T 2 + 6pa 2 + 3A 2 p 2 + 6 P A 2 + 3A 2 ) 
(p 2 a 2 +(i 2 + P A 2 + A 2 ) 3 / 2 

and excess excess cor — ^| cor — 3 = 



(25) 



, 2 (P 4 ^ 2 + 6p 2 a 2 + a 2 + 2A 2 p 3 + 6^V + 6pA 2 + 2A 2 ) ^ 
S ° (p 2 a 2 + a 2 + pA 2 + A 2 ) 2 
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Fig. 11. The product probability distribution with the 
correlation factor as a parameter, a) a 3D presentation; 
and b) sections for different correlation factors. 
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Fig. 12. The skewness (a) and the excess (b) for a product 
with 2 degrees of freedom and p as parameter. 



References 

Ardenne van, A. and Smits, F.M.A. 1996, in: Proc. of 
Workshop on Large Antennas in Radio Astronomy, 28-29 
February 1996, ESTEC, Noordwijk, The Netherlands, 59 

Barnbaum C, Bradley R. F. 1998, The Astronomical J., 115, 
2598 

Bregman J.D. 2000, Proc. of SPIE, 4015, 19 
Briggs F. H., Bell J. F., Kesteven M. J. 2000, Astrophysical J., 
in press. 

Cohen, J. 1999, Astronomy&Geophysics, 40, 6.8 
Deutsch, R. 1962, Nonlinear transformations of random pro- 
cesses, Prentice-Hall, Inc. Englewood Cliffs, N. J. , Ch. 4 
Ellington S. W., Bunton J. D., Bell J. F. 2000, Proc. of SPIE, 
4015, 400 

Fridman P. A., in: Proc. 8th IEEE Signal Processing Workshop 

on Statistical Signal and Array Processing, June 24-26, 

1996, Corfu, Greece, 264 
Fridman P.A. 1997, Radio Frequency Interference Rejection in 

Radio Astronomy Receivers, NFRA Note 664 
Fridman P.A. 1998, in: Proceedings of EUSIPC098 Ninth 

European Signal Processing Conference, Rhodes, Greece, 

8-11 September 1998, 2241 
Fridman P.A. 2000, in: Proceedings of the IEEE NORDIC 

Signal Processing Symposium, June 13-15, Kolmarden, 

Sweden, 375 

Kahlmann H.C. 1999, in: Review of Radio Science 1996-1999, 
ed. by W.Ross Stone, Published for the URSI by Oxford 
University Press, 751 

Lemmon J.J. 1997, Radio Science, 32, 525 

Leshem, A., van der Veen A. 1999 in: Perspectives on 
Radio Astronomy Technologies for Large Antenna Arrays, 
Smolders, A.B., and van Harlem, M.P., eds., ASTRON, 
Dwingeloo, 210 

Middleton D. 1972, IEEE Trans, on Electromagnetic 

Compatibility, EMC- 14 , 38 
Middleton D. 1977, IEEE Trans. on Electromaghetic 

Compatibility, EMC-19, 106 
Spoelstra T. A. Th. 1997, Tijdschrift van het Nederlands 

Elektronica- en Radiogenootschap, 62, 13 
Thompson, A. R. 1982, IEEE Trans, on Ant.&Propag., AP-30, 

450 

Whalen A. D. 1971, Detetection of Signals in Noise, Academic 
Press, N. Y. , ch. 4 

Weber, R., Faye, C, Biraud, F., Dansou, J. 1997, 
Astronomy&Astrophysics Suppl. Ser., 126, 161 

Widrow, B. and Stearns, S. 1985, Adaptive Signal Processing, 
Prentice-Hall, Inc., Englewood Cliffs, N. J. 

The Elizabeth and Frederick White Conference on Radio 
Frequency Interference Mitigation Strategies, December 
1999 CSIRO Radiophysics Laboratory Sydney Australia, 
|http: / / www.atnf.csiro.au /SKA /int mit /atnf / conf / 1 



When p approaches 1, both eqs. (25) and (26) reduce to 
eqs. (10) and (11), respectively. Fig. 12a, b show how the 
skewness and the excess depend on the RFI amplitude 
A with p as a variable. When p = 1 these curves are 
transformed into the curves in Fig. 2b, c corresponding to 
the sum of two squares of gaussian values. 
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